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We extend calculational techniques for static solitons to the case of field configurations with 
simple time dependence in order to consider quantum effects on the stability of Q-balls. These 
nontopological solitons exist classically for any fixed value of an unbroken global charge Q. We 
show that one-loop quantum effects can destabilize very small Q-balls. We show how the properties 
of the soliton are reflected in the associated scattering problem, and find that a good approximation 
to the full one-loop quantum energy of a Q-ball is given by uj — Eo, where ui is the frequency of the 
classical soliton's time dependence, and Eo is the energy of the lowest bound state in the associated 
' scattering problem. 
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INTRODUCTION 

A pure scalar theory in three dimensions with a cubic coupling can support classically stable, time-dependent 
non-topological solutions to the equations of motion carrying an unbroken global charge Q, called Q-balls [jjj. Super- 
^vq ' symmetric extensions of the standard model generically contain such objects They become particularly interesting 
in cosmological applications at small values of Q, because then it is easier for them to form in the early universe |$|. 
0^ ' In this regime, quantum corrections to the soliton's energy become increasingly important in determining its stabil- 
ity. The methods of jj] (see also earlier related work in [^) and theoretical justification in the Appendices of ||) 
provide an efficient, robust framework for computing quantum corrections to time-independent field configurations. 
In this Letter, we extend this approach to Q-balls. We show how to express the computation in terms of an effective 
y—{ ' time-independent problem. In this formalism, the full one-loop quantum correction can be computed efficiently. We 
also derive a very simple estimate for this result. The result is that we can compare the energy of the Q-ball in the 
quantum theory to the energy of free particles carrying the same charge Q, and determine if the Q-ball remains stable 
in the quantum theory. 

Our starting point is the classical analysis of Q-balls carried out in Jf|, ||. We will take the same simple model, 
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C=-(d ll <p){d»<p)-U(<p) (I) 



where ip is a complex field with unit charge under a global U(l) symmetry that is unbroken at ip = 0. We will consider 
<-h ' the potential 



U(<p) = ±M 2 \ V >\ 2 -A\<p\* + \\<p\ i . (2) 



A particular configuration (p(x, t) has charge 

Q = hJ d3x ^* dtip ~ vdt ^ • (3) 

Upon fixing the charge Q of the configuration via a Lagrange multiplier uj, we obtain the classical Q-ball solution 
as the minimum of the energy functional 

S w [ip] = f d 3 x^\d tV -iuj^\ 2 + I d 3 x Q|V^| 2 + C/ W (^ +ojQ (4) 



* e-mail: graham@physics.ucla.edu 



2 



with respect to independent variations of tp(x, t) and uj, where 
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U u {<p) = U(ifi) - -uj 2 cp 2 

The Q-ball solution then has simple time dependence 

<p(x,t) = e iwt ct>(x) 

so we must simply minimize 



(5) 
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by varying uj and (f>(x). As shown in Q|, if the quantity C/(0)/</> 2 is minimized at (f>o > 0, then for luq = y/2U(<po)/4 l o, 
the effective potential U Ua (<p) will have degenerate minima. For uj > u>q, a solution to the equations of motion is given 
by the bounce solution for tunneling in three Euclidean dimensions in the potential U u (0) . The bounce is the solution 
to 



with the boundary conditions 



lim (po{r) — and -r<Po( r ) 



= 0. 



(8) 



(9) 



r=0 



We can then find the solution using the shooting method detailed in JtJ, and minimize the resulting energy over to. For 
large enough Q, the optimal value of uj approaches ujq, allowing |jj to use the thin- wall approximation to demonstrate 
the existence of a global minimum, which is the Q-ball solution. For small Q, the optimal value of u> approaches 
M, and || uses the thick-wall approximation to show that there exists a global minimum in this case as well. Thus 
classically bound solitons exist all the way down to Q = 1. For all Q, we have ujq < uj < M at the minimum. We can 
thus consider the classical binding energy as a function of Q by comparing the Q-ball's energy to QM, the energy of 
a collection of free particles carrying charge Q. 

Although Q-balls are classically stable even as Q — > 1, the binding energy per charge is going to zero in this limit. 
This case is of particular interest for cosmological applications, however, since Q-balls of large charge, while favored 
energetically, are disfavored as the temperature increases by their low entropy. To answer the question of whether 
Q-balls have a significant chance of being formed in the early universe, we must therefore verify that the classical 
conclusions are not invalidated by quantum corrections. 



QUANTUM CORRECTIONS 



To compute the leading quantum correction the Q-ball energy, we extend the method of Q . We write the quantum 
field ip as the classical solution plus a quantum correction, which we write in corotating coordinates 

p(M) = e iwt (0 o (z)+r?(M)) (10) 

where we can then expand the quantum field r)(x, t) in small oscillations, which are given by the solutions to 

[(d t - iuj) 2 - V 2 + U"{ct> {x))} i>(x, t) = . (11) 

Parametrizing ip(x,t) = e luJt e~ lEt ip(x) gives the mode an energy E — uj, where the time-independent wavefunction 
ip(x) solves 



[-V 2 + U"(Mx))]i>(x)=E 2 i>(x) 



(12) 



which is an ordinary Schrodinger equation. The Casimir energy is given formally by the sum over zero-point energies 
of these oscillations 
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(13) 



3 



Since the spectrum of eq. ( [l2| ) is symmetric in E — > —E, we can sum over both signs of the energy and obtain 

£c"M~\ E {\E i +w\ + \E j -w\)= E max (M, 1^1) . (14) 



Ej>0 



We will use the methods of Q| to extract the quantum correction to the energy in terms of the continuum scattering 
data for the reduced problem of eq. (JT^). Since the potential is spherically symmetric, we can decompose the spectrum 
into partial waves I. We have wavefunctions 



Y tm (Q) 

ip e (x) = r)e(r) 

r 

for m = —£, — I + 1, . . . ,£ — 1, £. The radial wavefunction rjg(r) satisfies 

d 2 £(£+1) 
dr 2 r 2 



+ U"(Mr)) m{r)=E z m {r) 



(15) 



(16) 



with scattering boundary conditions. 

In each partial wave, we will find a continuum starting at E = M and possibly bound states with < Ej < M. 
(Since the spectrum is symmetric in E, we will only consider E > 0.) It is instructive to consider the properties of 
eq. ( |l2| ) that betray its origin from a field theory soliton. The full oscillation spectrum should have a zero mode in 
the £ = 1 channel, corresponding to the translation invariance of the Q-ball solution. The threefold degeneracy of this 
state corresponds to the three directions of translation. From eq. (|l^), we see that in the reduced problem, the zero 
mode appears as a bound state with energy E — to. Since this state appears in the £ = 1 channel, there must exist 
an even more tightly bound state in the I = channel. In the case of an ordinary static solution, this state would 
correspond to an instability of the full soliton. But from eq. ([l4]), we see that the destabilizing effect of this mode is 
neutralized by the time dependence of the classical solution, which results in it making the same contribution to the 
mode sum as the zero modes do. All other modes have energies greater than to. 

Having rewritten the Casimir energy in terms of the eigenmodes of the reduced scattering problem in eqs. ( |l2| ) and 
(p"6|), we are prepared to apply the methods of [Q. We obtain the renormalized Casimir energy as a sum over partial 
waves t. In each partial wave we have a sum over positive energy bound states Ej t £ and an integral over continuum 
states weighted by the density of states p(k), where k = \/E 2 — M 2 . We subtract the corresponding integral in the 
free case, obtaining 



(=0 



f°° dh 

/ — E ( Pe (k) - P ° e (k)) +Vmax( Wl %) 



(17) 



This integral diverges in the unrenormalized theory, because we have not yet included the contribution of the countert- 
erms. We will compute the continuum integral by relating the difference between the interacting and free densitites 
of states to the phase shifts, 



pt {k)-p%k) = ^5 t {k) 



(18) 



The Born expansion for the phase shift is then in exact correspondence with the expansion of the effective energy in 
one-loop diagrams with all possible insertions of the background field. Subtracting the first two Born approximations 
from the phase shift corresponds to subtracting the first two diagrams in this expansion, which are the only divergent 
terms. The remaining integral is then finite. We then add the divergent terms back in, together with counterterms, 
as ordinary Feynman diagrams. Full details are given in pi. 
We thus obtain the renormalized Casimir energy 



£ c [<M = r(% ] + ]T(2^ + i) 



e=o 



i°° v E i ( Se{k) ~ ^ ~ s " 2) (fc) ) + ^ max(cj ' Ej 



(19) 



6 { />(k) and 5Y>{k) are its first and second Born approxi- 



where 8i(k) is the scattering phase shift in partial wave P '" ! A '~ ; 

mations, and r' 2 ' [<j>o\ is the contribution to the energy from the two-point function, computed in ordinary Feynman 
perturbation theory. This piece includes the counterterms, which we fix using physical renormalization conditions. 
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We demand that the tadpole graph vanish, and that the mass of the free <p particle is unchanged. In addition to 
holding the location of the pole in the propagator fixed at M, we also perform wavefunction renormalization so that 
its residue is unchanged as well. 

We can now compute eq. (19) directly. We simply require the scattering phase shifts, their Born approximations, 
and the bound states of eq. (12). Efficient algorithms for obtaining these are detailed in Q. The contribution to the 
energy from the two-point function is computed using conventional techniques, giving 



4q 2 dq 
(4tt)* 



arctanb 



5tt 



V? 2 + 4M 2 3^3 



+ i \x(q)\ 2 -^ 



( 2tt 



(20) 



where x(q) and <j)(q) are the spatial Fourier transforms of U"(<f>o(r)) — M 2 and <f>(f) respectively. 

Examining this calculation in detail yields a very accurate estimate for the quantum correction to the energy, which 
is very easy to compute. Using the analysis of bound states above, we can separate eq. ( JTgj ) into 



£c[<h]= T 



(2)r 



5>- 

e=o 



+ (W - E ) 



(21) 



where Eq is energy of the most tightly bound state, which appears in the I — channel. It is the only state with 
energy less than lo. We define the reduced Casimir energy as 



£c"[<M = r(% ] + ^(2^ + i) 

1=0 
oo 

= r<% ] + E(2* + i) 



dk d 
7r dk 



dk 
o n 



(5 t {k)-5f\k)~8f\k))+Y^E } 

3 

(E M)^ (S e (k) - S^(k) - 5f >(*)) + - M) 



(22) 



where we have used Levinson's theorem in the second line. This quantity is simply the Casimir energy of a time- 
independent soliton giving rise to the reduced small oscillations of eq. (|l^). (Of course, such a soliton would not solve 
the field theory equations of motion, but we could imagine holding it in place with an external source). The reduced 
potential is shallow and slowly varying, especially in the limit of small Q, which corresponds to lo approaching M. 
It causes only a slight deformation of the small oscillations spectrum — in particular, there is only one state bound 
more tightly than lo. For a generic potential of this kind, the contributions from the bound states and continuum will 
be opposite in sign; roughly, rearrangement of the continuum spectrum partially compensates for the effect of the 
states that become bound. 

A direct application of j| allows us to evaluate the full result of eq. (22). We can also estimate this result in the 
derivative expansion. To lowest order, we have simply the effective potential contribution 



Crcduccd.DE r / 1 

' l<Po\ 



where 



^M 4 ((l + z) 2 log(l + ,)-,-^ 2 



U"(Mr))~M 2 



M 2 



(23) 



(24) 



Using either technique, explicit computations show that this reduced Casimir energy is very small compared to the 
classical binding energy of the Q-ball (typically 5% or less for small Q). Thus we lose very little accuracy by dropping 
this term, obtaining a very simple estimate for the Casimir energy: 



£c[(M 



&q] = U - Eq . 



(25) 



APPLICATIONS 



To see whether the Q-ball is stable, we must compare its energy to the energy of a state with the same charge built 
on the trivial vacuum 



B[(po] = £[0o] - QM 



(26) 
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FIG. 1: Q-ball binding as a function of Q, in units of M. Parameters are A = 0.325M and X = 0.055 (left panel), and 
A = 0.425M and X — 0.095 (right panel). Shown are three calculations of the difference B between the energy of a Q-ball 
and the energy of a state with charge Q built on the trivial vacuum: the classical approximation B dlMB [<j>o], the full one-loop 
calculation Bmi[4>q], and the estimated one-loop result B cst [4>o]- 



Figure [j] shows the result of different calculations of £ , each as a function of Q, for two choices of the coupling 
constants. In both cases, they are chosen so that ip = remains the global minimum of U(ip). 1 We work in units of 
M, which sets the scale of the problem. In the classical approximation, 

Bc„W=^N-QM (27) 

we see the result of ||: the Q-ball is stable for all Q, though the binding energy per charge is going to zero as Q — > 0. 
In the full one-loop calculation, 

S M [<fo] = £ u [<fo] + £c [<fo] ~ QM (28) 

we see that the quantum corrections overwhelm the weak classical binding up to Q min « 7. Above this value, the 
Q-ball is stable. Finally, we see that using eq. ([25| ) to approximate to the one-loop result by taking 

£L t [4> ] - e u [M + &t [<M - QM (29) 

yields a result that is very close the full one- loop result. It is also interesting to note that this approximation is 
particularly good near the value of Q at which the Q-ball becomes bound in the full one-loop calculation. 



CONCLUSIONS 



We have computed one-loop quantum corrections to the energies of Q-balls. For small Q, these corrections can play 
an important role in determining the stability of these extended objects. Since Q-balls are genuine solutions to the 
classical equations of motion, and we are working in a regime where the coupling constants are small, the one-loop 
approximation should be very good. We have seen that the one-loop correction is indeed very small compared to the 
classical energy of the Q-ball, although can be significant when compared to its binding energy, the difference between 
the energy of the Q-ball and the energy of a collection of free particles carrying the same charge. The higher-loop 
corrections should be correspondingly smaller than the one-loop corrections, and therefore can safely be neglected. 
Of course, if the theory contains other particles coupled to the ip field, we may need to include their contributions as 
well, using the same techniques as we have developed here. For small coupling constants, the one-loop approximation 
should continue to be reliable. Within the one-loop approximation, we have seen that the quantum correction can be 
very accurately estimated at small Q by considering the difference between the Lagrange multiplier ui, which gives 



1 I thank M. Postma for reminding me of this requirement. 



6 



the frequency of the Q-ball's time dependence, and the energy Eq of the lowest bound state of the small oscillations 
potential U"((j>o(x)). The result is a prediction of the minimum value of Q for which the Q-ball is stable in the 
quantum theory. For typical values of the coupling constants, we find Q min 7. 
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